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1. Introduction 

Carbon nanotubes show great promise for use as nanoscale materials due to their 
high tensile strength in the axial direction and remarkable flexibility in bending pQ. 
Another important characteristic of carbon nanotubes is their high flexibility in the 
radial direction. In fact, the magnitude of radial stiffness of an isolated carbon nanotube 
is considerably less than that of axial stiffness [2], which allows a reversible change 
in the cross-sectional shape on applying a hydrostatic pressure. Such a pressure- 
induced radial deformation and the associated change in vibrational modes are useful 
to probe the structural properties of nanotubes. More interestingly, this deformation 
alters other physical properties of nanotubes and related materials, such as electronic 
[HIHE1E1H1E1E] and optical pU QU M, DSl M, IS] properties. For instance, application 
of hydrostatic pressures have been found to induce drastic changes in the electrical 
conductance of nanotubes [161 E]j implying the relevance of radial deformations in 
carbon nanotube applications. 

Thus far, many experimental and theoretical studies have been carried out on radial 
deformations of carbon nanotubes induced by hydrostatic pressures [HI HH EE} 
EaESESEElEniESE^ most of which have focused on single- 

walled nanotubes (SWNTs) and their bundles. Successive investigations have revealed 
flattening and polygonalization in the cross section of SWNTs under pressures of the 
order a few GPa [TOl [19], in which the stable cross-sectional shape is determined by 
minimizing the elastic energy of the cylindrical tube under constraints. Employing 
these transformation properties for developing nanoscale pressure sensors has also been 
suggested [3U [35] . 

As compared to the intensive studies carried on SWNTs, studies on the radial elastic 
deformation of multiwalled nanotubes (MWNTs) are lagging behind. An important 
feature of MWNTs is that they consist of a set of concentric graphite cylindrical walls 
mutually interacting via the van der Waals (vdW) forces. This core-shell structure 
produces an encapsulation effect under external pressure, wherein the outer walls 
collectively function as a protective shield. In principle, the encapsulation effect 
enhances the radial stiffness of MWNTs. However, this is less obvious if the number 
of concentric walls, N, is much greater than unity. For the latter case, outside walls 
possess large diameters so that the application of a hydrostatic pressure leads to a 
mechanical instability in the outer walls. This instability is, nonetheless, compensated 
by the relative rigidity of the inner walls with small diameters. These two competing 
effects imply the possibility of a new cross-sectional shape transition of MWNTs induced 
by hydrostatic pressure, whereas such the transition is still to be explored. 

In this article, we demonstrate a novel radial deformation, called the radial 
corrugation, of MWNTs with N ^> 1 under hydrostatic pressure. In the corrugation 
mode, outside walls show significant radial corrugation along the circumference, while 
the innermost tube maintains its cylindrical symmetry (Fig. Hfb)). We found various 
corrugation modes can be obtained by tuning the innermost tube diameter D and the 
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number of constituent walls N, which is a direct consequence of the core-shell structure 
of MWNTs. These results provide useful information for developing nanoffuidic 
[36| 1371 I3~8| 139] or nanoelectrochemical [UH E] devices whose performance depends 
crucially on the geometry of the inner hollow cavity of nanotubes. 

2. Method 

2.1. Mechanical energy of MWNT 

The stable cross-sectional shape of a MWNT under a hydrostatic pressure p can be 
evaluated by using the continuum elastic theory for cylindrical shells. The effectiveness 
of the continuum approximation for modeling MWNTs has been demonstrated in a 
series of studies previously [121 H31 HH S3 HE], particularly for multiwalled nanotubes 
containing a large number of carbon atoms. In the continuum approximation, the 
mechanical energy U of a MWNT per unit length in the axial direction is the sum of the 
deformation energy Ud of all concentric walls, the interaction energy Uj of all adjacent 
pairs of walls, and the potential energy Q of the applied pressure. All the three energy 
terms are functions of p and the deformation amplitudes Ui(p,6) and Vi(p,9), which 
denote the radial and circumferential displacements, respectively, of a surface element 
of the zth wall (see Fig. [2]). Thus, U can be written as 



It should be noted that our theoretical model assumes a very long, straight MWNT 
with both ends being free. On applying hydrostatic pressure, therefore, the tube is 
freely (but slightly) elongated in the longitudinal direction because of the absence of 
reactive force at the ends. This means the absence of longitudinal deformation in our 
theoretical condition, which allows us to consider only the cross sectional deformation 
under hydrostatic pressure. 



U = U[p, mip, 0),Vi(p, 9)] =U D + U I + Q. 



(1) 




Figure 1. (color online) Cross-sectional views of (a) elliptic (n = 2) and (b) corrugated 
(n = 5) deformation modes observed for iV-walled nanotubes with (a) N — 29 and 
(b) N — 30. The innermost tube diameter D = 5.0 nm is fixed. The mode index n 
indicates the wave number of the deformation mode along the circumference. 
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Figure 2. (color online) Illustrations of geometric parameters of our continuum elastic 
model. Left: Sketch of cross section of MWNT subjected to hydrostatic pressure 
p. Here, indicates the radius of the ith cylindrical wall of the MWNT. Right: 
Sketch of displacement of surface element of ith wall. Here, u(p, 0) and v(p, 9) denote 
the deformation amplitudes of the surface element in the radial and circumferential 
directions, respectively. 



An explicit form of the function (CQ) is obtained as follows. Firstly, Ud of all 
concentric walls is expressed as 



N 

•i=i 
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Ui + Vi H 
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+ f Ut y }<">■ (2) 



where u' = dujdQ. In Eq. ([2]), the first term in the curly brackets describes the stretching 
energy of the walls in the tangential direction and the second term describes the bending 
energy (The derivation of the formula (T5]) involves some complex calculations and is 
presented in Appendix.) The two parameters cti and /3j characterise the mechanical 
stiffness of the walls for in-plane stretching and bending, respectively They are defined 
as follows: 

Eh , „ Eh 3 

ati = j and A = — — — , 

1 — v l 12 (1 — v 1 ) 

where E is Young's modulus (E = 1 TPa), v is Poisson's ratio [v = 0.27), and h 
indicates the thickness of individual walls (h = 0.34 nm). Hereafter, r, indicates the 
radius of the ith wall in the absence of pressure; thus D = 2r\. The spacing between 
adjacent walls is set to be |r, — rj ±1 | = 0.344 + 0.1e~ D / 2 nm according to Ref. [4T] . 
Secondly, the explicit form of Ui is given as 

u i = J2 ''' ' / (ui-u i+1 ) d9 + J2 l '\ / (ut-Ut-i) dO, (3) 
T~\ 2 Jo ~o 2 Jo 
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where the vdW interaction coefficients are functions of rj and rj and given as follows 
03: 

AOOIttcct 12 mOvrecx 6 ^ 



Here we set 



n / 2 dO , 4r,r 



FT? = — / and AT, 



in which a denotes the chemical bond length between neignbouring carbon atoms within 
a layer (a = 0.142 nm), and e and a are the parameters that determine the vdW 
interaction between two layers (e = 2.968 meV and a = 0.3407 nm) [48J. 

Finally, we consider an explicit form of Q. Since a cylindrical shell subjected to p 
is a conservative system, the change in the potential energy as the cross section of the 
MWNT deforms is the negative of the work done by the pressure during the deformation. 
Hence, we have 

n = -p(nr* N -S*), (4) 

where S* is the cross-sectional area after deformation. A simple calculation yields the 
following expression: 

n [ 2n ( , u% + v n ~ Un'vn + u N v N ' \ 

11 = p \r N u N -\ \dd. (5) 



o \ 2 
Refer Appendix for the derivation of the formula (JSj). 

2.2. Critical pressure p c and deformation mode n 

Our aim is to evaluate the critical pressure p c above which the circular cross section 
of MWNTs is elastically deformed into non-circular one. To carry out the analysis, we 
decompose the radial displacement terms as Ui (p,0) = u f\p) +5ui(9). Here, uf\p) 
indicates a uniform radial contraction of the zth wall at p < p c , whose magnitude is 
proportional to p. The Sui(9) describes a deformed, non-circular cross section observed 
just above p c . Similarly, we can write Vi(p, 6) = 5vi(6), since v[ (p) = at p < p c . 

Applying the variation method to U with respect to Uj and we obtain a system 
of 2N linear differential equations given by 

ai(5ui + dvi - 7»^) + PiT)'" + p5 i:N (5ui + 5v-) 
+ (1 - S iiN ) c iii+1 ri (Sui - 5u i+ i) 

+ (1 - 6i,i) c^i-tri (Sui - Sui-i) = 0, (i = 1, • • -, N) (6) 

and 

ai(5ui+5vi , ' + j i rk)-p i 7i?+p5 itN (5u i '-5v i ) = 0, (i = l,---,N) (7) 

where ji = uf\p)/ri and rji = u^—Vi. In deriving Eqs. (jSJ) and (J7J) through the calculus of 
variation, quadratic or cubic terms in 5ui and 5vi were omitted since we consider elastic 
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deformation with sufficiently small displacements. In addition, the terms consisting only 
of uf^ and p are also omitted; the sum of such terms should be equal to zero since 
represents an equilibrium circular cross-section under p\^ 

Clearly, 5ui, 5vi and their derivatives are periodic in 9. Hence, general solutions of 
Eqs. (E]) and (J7J) are given by the Fourier series expansions 

oo oo 

5uj(9) = Sp,j(n) cosnO and Svj(9) = ^ Suj(n) sinnfl. 

n=l n=l 

Introduction into Eqs. ([6]) and (J7|) results in the matrix equation Cu = 0, in which 
the vector u consists of 5pi(n) and 8ui(n) with all possible i and n, and the matrix C 
involves one variable p as well as parameters such as Oj, /3j, • • • etc. It should be noted 
that, due to the orthogonality of cos n6 and smn9, the matrix C can be expressed by a 
block diagonal matrix of the form C = C n= \ © C n=2 © • • •• Here, C n=m is a 2N x 27V" 
submatrix that satisfies C n=m u n=m = 0, where u n=m is a 2iV-column vector composed 
of Sp,i(n = m) and 5i>i{n = m). As a result, the secular equation det(C) = that 
provides nontrivial solutions of Eqs. and ([7]) is rewritten by 

det (C n= i) det (C n=2 ) • • • = 0. (8) 

Solving Eq. (jSJ) with respect to p, we obtain a sequence of discrete values of p each 
of which is the smallest solution of det (C n=m ) ■ ■ • = (m = 1, 2, ■ ■ ■). Among these 
values of ps, the minimum one serves as the critical pressure p c that is associated with 
a specific integer n = m. From the definition, the p c associated with a specific m 
allows only 5fii(n = m) and bi>{{n = m) be finite, but it requires Sfli(n ^ m) = and 
bviin ^ m) = 0. Immediately above p c , therefore, the circular cross section of MWNTs 
becomes radially deformed as described by 

Ui{6) = u\ ° (p c ) + #A*i( n ) cos 72$ and Vi(6) = 5ui(n) sin n6, 

where the value of n is uniquely determined by the one-to-one relation between n and 

Pc- 



3. Radial corrugations of MWNTs 

3.1. Cross-sectional view 

Figure Q] illustrates cross-sectional view of two typical deformation modes — (a) elliptic 
(n = 2) and (b) corrugation (n = 5) modes — of a MWNT with D = 5.0 nm. For 
Fig. [Ua) and 1(b), iV = 29 and 30, respectively. In the elliptic mode, all constituent 
walls are radially deformed. On the contrary, in the corrugation mode, outside walls 
exhibit significant deformation, while the innermost wall maintains its circular shape. 
The deformation mode observed just above p c depends on the values of iV and D for 
the MWNT under consideration. It will be shown below that larger iV and smaller D 
favor a corrugation mode with larger n. 

| The fact that the sum equals zero determines the function form of (p) . 
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Figure 3. (color online) (a) Critical pressure curves showing p c required to produce 
radial deformation of N- walled nanotubes with fixed D. The mode index n that 
characterizes the radial deformation mode observed just above p c is also shown, (b) 
Branches of solutions p(N) for secular equation det(C) = (refer text). The iV- 
dependence of p(N) for each deformation mode n is displayed. The innermost tube 
diameter is set to be D = 4.0 nm for all curves. For a fixed N, the minimum value of 
p among the branches functions as the critical pressure p c (N) just above which radial 
deformation takes place. The vertical arrows indicate the phase boundary points across 
which the deformation modes change. 



3.2. Critical pressure curve 

Figure [3] (a) shows J) c clS 8b function of iV for various values of D. The mode index n of 
the deformation mode observed just above p c for fixed N and D is also shown. For all 
D, p c increase with iV followed by a slow decay, except for the case of D — 8.0 nm. The 
increase in p c in the region of small iV is attributed to the enhancement of radial stiffness 
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of the entire MWNT by encapsulation. This stiffening effect disappears with further 
increase in N, resulting in the decay of p c (N). A decay in p c implies that a relatively low 
pressure becomes sufficient to produce radial deformation, thus indicating an effective 
"softening" of the MWNT. Such a decay is also observed for D = 8.0 nm and larger 
D, in principle, if a sufficiently large N is considered (but omitted in Fig. El (a)). The 
two contrasting behaviors, stiffening and softening, are different manifestations of the 
encapsulation effect of MWNTs. 

We observe in Fig. [3] (a), only the elliptic deformation is realised in the stiffened 
region, while certain radial corrugations with n > 3 are formed in the softened region. 
The observation of the former phenomenon can be understood by considering the radial 
deformation of a single-walled cylindrical shell. In the stiffened region, since iV is very 
small, the relation \r± — rjy| <C rjy holds. This implies that in this region, a MWNT 
behaves effectively as a single-walled cylindrical shell with wall thickness \r\ — r^\. For 
a single- walled shell, the secular equation det(C) = gives solutions of 

Eh 3 n 2 - 1 

P ~ 12(1-^)7^ l + (h?/12r%) ^- 2 > 3 '"--)- 

Hence, p always has the minimum value when n — 2. This scenario applies to MWNTs 
if N is sufficiently small, resulting in the elliptic mode of n = 2 in the stiffened region. 

The most striking observation is the successive transformation of the cross section 
with an increase in N. For D = 4.0 nm, for example, the deformation mode observed 
just above p c jumps abruptly from n = 2 to 4 at N = 17, followed by successive 
emergences of higher corrugation modes with larger n. The critical number of walls 
N = N c separating the elliptic phase from the corrugation phase is identified to iV 
that yields a cusp in the curve of p c (N). In contrast, no singularity is observed in 
the curve of p c (N) at values of N, which separate two neighboring corrugation phases. 
We emphasize that at these phase boundaries, one additional wall induces a drastic 
change in the cross-sectional shape of the MWNT under the hydrostatic pressure. These 
transitions in n originate from the two competing effects inherent in MWNTs with 
N ^> 1, that is, the relative rigidity of the inner walls and the mechanical instability of 
outer walls. A large discrepancy in the radial stiffness of the inner and outer walls gives 
rise to a maldistribution of the deformation amplitudes of concentric walls interacting 
through the vdW forces, which consequently produces an abrupt change in the observed 
deformation mode at some N. 

Figure [3](b) explains the absence of the corrugation mode of n = 3 in MWNTs 
with D = 4.0 nm. This figure shows the iV-dependence of the solutions p(N) for the 
secular equation det(C) = 0. As mentioned earlier, the secular equation gives various 
values of p, each of which is associated with a specific mode index n. Among the values 
of p, the minimum value gives the critical pressure p c just above which cross-sectional 
deformation takes place. Figure [3]^b) depicts the ^-dependence of p(N) for several n 
values, where the innermost tube radius is fixed to be D = 4.0 nm. For iV < 17, the 
values of p for n = 2 are less than those for n > 3, which implies that the elliptic mode 
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Innermost tube diameter: D [nm] 

Figure 4. (color online) Phase diagram of radial deformation modes observed above 
p c . Various corrugation modes are obtained depending on the values of N and D. 

occurs for MWNTs with iV < 17. However, for N > 17, the minimum p corresponds 
to n = 4, implying the occurrence of the corrugation mode of n = 4. It should be 
noted that for n = 3, p can never attain the minimum value at any N. This is why 
the corrugation mode of n = 3 cannot be observed for MWNTs with D = 4.0 nm. A 
parallel discussion accounts for the absence of the modes of n = 3, 4 for D = 5.0 nm 
and the corrugation modes for D = 8.0 nm within the range of N we have considered. 
It also follows from Fig. [3](b) that the cusps in the curves p c (N) occur only at the phase 
boundary N c separating the elliptic phase (n = 2) from a corrugation phase (n > 3), 
while no singularity appears at the boundaries of N between neighbouring corrugation 
phases. 

3.3. Phase diagram 

Figure H] shows a phase diagram of the radial deformation modes in MWNTs observed 
above p c . The bottom region below a chain of solid circles (colored in gray online) 
corresponds to the elliptic phase, and the top regions surrounded by other symbols 
indicate the corrugation phases associated with n. It is clearly observed that above 
the elliptic phase, multiple corrugation modes are formed depending on the values of 
iV and D. It is also observed that smaller D and larger iV favor corrugation modes to 
the elliptic modes. Furthermore, larger iV yield higher corrugation modes with larger n. 
To our knowledge, this is the first demonstration of the phase diagram of the multiple 
corrugation transitions, while conditions for the higher modes to occur were previously 
argued [23]. A quantitative examination of the present results by sophisticated atomistic 
simulations [321 EO, El] may provide a better understanding of the corrugation properties 
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10 20 30 40 50 1 5 10 50 
Index of walls: Index of walls: / 



Figure 5. (color online) Left: Deformation amplitude of each zth concentric wall of 
MWNT with N — 50. The mode index n of the associated deformation mode is also 
indicated. Right: Double-logarithmic plots of data plotted in left-hand-side figure. 



of MWNTs. 

4. Persistent cylindrical geometry of the innermost tube 

We have mentioned that in all corrugation modes, the innermost tube maintains its 
cylindrical symmetry. This persistence of cylindrical symmetry can be demonstrated by 
plotting the deformation amplitudes Sfli. Figure [5] shows the normalized deformation 
amplitudes, |5/2j/5/2/v|, of individual concentric walls for the MWNT with iV = 50. The 
innermost and outermost walls correspond to i — 1 and % = 50, respectively, and the 
associated n is also shown. Double-logarithmic plots of the same data for the corrugation 
mode (n = 6) are also presented in the right-hand-side figure. The figures clearly show 
that in the corrugation mode, the deformation amplitudes of the innermost tube are 
significantly smaller than those of the outer walls by several orders of magnitude. The 
persistence of cylindrical symmetry of the innermost tube will be useful in developing 
nanotube-based nanofluidic [361 E3 [381 EH] or nanoelectrochemical devices 001111], since 
both utilize the hollow cavity within the innermost tube. In fact, several different 
types of intercalated molecules such as diatomic gas, water, organic, and transition 
metal molecules are known to fill the innermost hollow cavities of nanotubes [37] and 
exhibit various intriguing behaviors that are distinct from those of the corresponding 
bulk systems [521 [531 El]- For the intercalates, the innermost tube of MWNTs with 
N ^> 1 functions as an ideal protective shield, since it maintains its cylindrical geometry 
even under high external pressures up to an order of GPa. 
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5. Discussions 

We have observed that radial corrugation modes yield a large difference in the 
deformation amplitudes between the outer and inner walls. Hence, the modes disturb 
the equal spacings between the concentric walls in MWNTs. This possibly affects 
electronic and vibrational properties of the entire nanotube, thus triggering a change in 
its electronic and thermal conductance. Such pressure-induced changes, if they occur, 
are of practical use for developing MWNT-based pressure sensors [3U [35]. Besides, 
these changes should be relevant to the characteristics of carbon nanotube composites 
[551 [56] that contain many clusters of MWNTs dispersed in an elastic medium. In the 
composites, the vdW interaction between adjacent tubes is sufficiently strong to yield 
radial corrugation, which affects the mechanical strength or the conducting properties 
of the composites. Intensive studies on these issues are expected to yield new methods 
for engineering next-generation devices and nanomaterials. 

Before presenting the conclusion, we mention the relevance of the lattice registry 
in radial corrugation phenomena. The atomic lattice registry (i.e., the degree of 
commensurance in atomic structures between neighbouring carbon layers) is known to 
play a prominent role in determining the optimal morphology of fully collapsed MWNTs 
[571 EH [59]. For a similar reason, it possibly affects the cross-sectional shape of the 
corrugation modes, if the interlayer spacings partially vanish by applying pressure much 
higher than that we have considered; Crumpling or twisting in outside walls may be 
observed depending on the atomic configuration of the MWNT. We also conjecture 
that the difference in the lattice registry between adjacent layers will induce a shift 
in the phase boundary depicted in the phase diagram of Fig. [H For instance, the 
AA stacking of carbon atoms [59j in the radial direction will lead to an effective load 
transfer between neighbouring layers. Thereby, rigid inner walls can support effectively 
the mechanical instability of outer walls. This implies that the phase boundary curve 
shown in Fig. H] separating the elliptic phase (n = 2) from corrugation phases (n > 3) 
shifts downward. On the other hand, the AB stacking tends to shift upward toward the 
phase boundary since the rigid inner walls provide minimal support. The magnitude of 
the phase boundary shift will be quantitatively determined in our future studies. 

6. Conclusion 

In conclusion, we have demonstrated the presence of multiple radial corrugations 
peculiar to MWNTs under hydrostatic pressures. Theoretical investigations based on 
the continuum elastic theory have revealed that MWNTs consisting of a large number 
of concentric walls undergo elastic deformations at critical pressure p c ~ 1 GPa, above 
which the cross-sectional circular shape becomes radially corrugated. A phase diagram 
has been established to obtain the requisite values of N and D for observing a desired 
corrugation mode. It is remarkable that in all corrugation modes, the cylindrical 
symmetry of the innermost tube is maintained even under high external pressures. This 
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Figure 6. (color online) Part of cross section of ith cylindrical wall undergoing radial 
deformation. Cross-sectional deformation leads to the translation of the shaded region 
colored in dark green (located at the point fa, 9)) into the adjacent shaded region 
colored in light green. The vectors n and n' indicate the normals to the centroidal 
surface (dashed-dotted) before and after deformation, respectively. The shaded region 
after deformation exhibits stretching in the circumferential direction with rotation by 



persistence of the cylindrical symmetry of the innermost tube of MWNTs is completely 
in contrast to the pressure-induced collapse of SWNTs. We believe that a study of this 
behavior of MWNTs will shed light on the potential of MWNT-based devices. 
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Appendix 

Appendix A.l. Derivation of the deformation energy Up 

In this section, we provide the explicit derivation of Up introduced in Eq. (TSJ). Consider 
the ith cylindrical wall of a long and thin circular tube with thickness h. A surface 
element of the cross-sectional area of the wall is expressed by (rid8)dz, where 9 is the 
circumferential angle around the cylindrical axis and z is a radial coordinate measured 
from the centroidal surface (indicated by A in Fig. [6]). The stiffness of the element for 
stretching along the circumferential direction is given by E/(l — z/ 2 ), where E and v are 
Young's modulus and Poisson's ratio, respectively, of the wall. Thus, the deformation 
energy U$ of the ith wall per unit length in the axial direction is written as 




(A.l) 
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Here, e(z, 9) is the extensional strain of the element in the circumferential direction. For 



the calculation, we adopt the strain-displacement relation introduced by Sanders [60] 
subject to the the following two assumptions: (i) strains are small and (ii) the rotation 
angles (3 of the normal to the centroidal surface are 'small but finite'. Then, we obtain 
the relation 1 60 1 



s(z,9) = e(9) + zk(9), (A.2) 
where e and k are expressed in terms of Ui(9) and Vi(6) as follows: 

Ui + v / 1 / u'; — vi \ 2 , uJ' — vJ , . „ s 

^ = ^r^ + o(^r^) and k = -^t^- ( A -3) 



The relation (1A.2j) states that the circumferential strain e of a volume element can be 
decomposed into in-plane stretching e{9) and bending-induced stretching zn{9). The 
presence of the latter is implicitly illustrated in Fig. [BJ where the outermost surface is 
more elongated in the circumferential direction than the innermost surface due to the 
thickness of the wall. 

From (1A.1I) and (1A.2j) . we obtain 



Ehr- r 2n F!h 3 r- r 27r 

(Y) _ rjfll i l 2 Hill ft / 2 



= / e *d0 + LL- / k d9. (A.4) 

D 2(1- i/ ) Jo 24 (1 - 2/ 2 ) 7o V ; 



By substituting ( 1A.3I) into ( 1A.4j) . followed by some calculations, we obtain the desired 
form of Up — J2iLi as given in (j2J). 

Appendix A.2. Derivation of the potential energy Q 

This section presents the derivation of the potential energy Q of applied pressure given 
in Eq. (jSJ). Suppose that an elastic cylindrical shell having a radius r N is subjected to a 
hydrostatic pressure p. Since fl is the negative of the work done by the external pressure 
during cross-sectional deformation, it is expressed as 

Q = -p(rcr 2 N -S*), (A.5) 

where S* is the cross sectional area of the shell after the deformation (the sign of p is 
assumed to be positive inward). 

We now consider the explicit functional form of S*. We denote the Cartesian 
coordinates of a point on the circumference of the cross section by (x*,y*). After the 
deformation, they are parametrised by the angle 9 as follows: 

x*{9) = [rjy + un(9)} cos 9 — vn{9) sin 9, 

y*(9) = [r N + u N (9)} sin 9 + v N (9) cos 9. (A.6) 

Here, un{9) and v^(9) are the deformation amplitudes of a surface element of the shell in 
the radial and circumference directions, respectively. S* is obtained by the line integral 
around the circumference C as follows: 

S* = ~<f c (-y*dx* + x*dy*), 
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or equivalently from, 




d9 d9 
Substituting KKM into (TA~7j) yields 



v ~nr + x -nr\ de - ( A - 7 ) 



S* = irr 2 N + ^J o \%r N u N + u 2 N + v 2 N - ~^- v n + u N~j]f ] d ®> ( Ai 
where the periodicity relation 



f 2n dv N 



o 



d e de = 

is employed. From (1A.5I) and ( 1A.8I) . we finally obtain the desired result 

n f 2n ( t u n + v n~ u n'v n + u N v N '\ 

Jo \ TnUn 2 J ^ ' 

where u' = du/de. 
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